数値微分は、微積分の無限に滑らかな性質から、離散的で有限なデジタル計算の世界への重要な転換を示します。無限小の極限ではなく、測定可能なステップサイズ $h$ を用います。関数 $f$ の $x_0$ における理論的な導関数は $$f'(x_0) = \lim_{h \to 0} \frac{f(x_0 + h) - f(x_0)}{h}$$ として定義されますが、コンピュータシステムは直接極限を計算できません。代わりに、有限差分公式を使用し、明確に定量化できるペナルティである 切り捨て誤差。
1. 導関数の幾何学的性質
$f'(x_0)$ を近似するためには、隣接する点に注目します。方向の選択によって、2つの主要な公式が導かれます：
- 前進差分公式： $h > 0$ の場合に使用します。$x_0 + h$ の先を見るように設計されています。
- 後退差分公式： $h < 0$ の場合に使用します。$x_0 + h$（ここで $h$ は負）を逆に見るように設計されています。
実際の工学問題、たとえば曲線の軌道の弧長を計算する場合など、これらの近似に頼ることが多いです：$$L = \int_{0}^{48} \sqrt{1 + (f'(x))^2} dx = \int_{0}^{48} \sqrt{1 + (\cos x)^2} dx$$ 関数 $f(x)$ が離散的なセンサー点でのみ分かっている場合、数値微分が唯一の解決策となります。
2. 求積法による数学的導出
$f'(x_0)$ を近似するために、まず $x_0 \in (a, b)$ かつ $f \in C^2[a, b]$、$x_1 = x_0 + h$ であると仮定します。$x_0$ と $x_1$ を通る最初のラグランジュ多項式 $P_{0,1}(x)$ を構築します：
ステップ1：補間関数の構築
$f(x) = P_{0,1}(x) + \frac{(x - x_0)(x - x_1)}{2!} f''(\xi(x))$
ステップ2：微分
両辺を微分し、$x = x_0$ で評価することで基本的な関係が得られます：
$$f'(x_0) = \frac{f(x_0 + h) - f(x_0)}{h} - \frac{h}{2} f''(\xi)$$
3. 誤差項と収束性
項 $-\frac{h}{2} f''(\xi)$ が私たちの切り捨て誤差です。この公式は精度が $O(h)$であることを示しています。つまり、ステップサイズ $h$ を半分にすると、誤差も概ね半分になります。しかし注意が必要です：$h$ が小さくなるほど切り捨て誤差は減少しますが、最終的には分子におけるほぼ等しい数値の差引きにより、 丸め誤差 分子におけるほぼ同じ数値の減算によって、増加します。
🎯 核心原則：有限差分
数値微分は極限を有限の弦に置き換えます。我々の近似の品質は、ステップサイズ $h$ と関数の滑らかさ（2次導関数）に厳密に依存します。
$f'(x_0) \approx \frac{f(x_0+h) - f(x_0)}{h}$ であり、誤差の上限は $\frac{h}{2} \max|f''(\xi)|$ です。